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This paper is motivated by the suggestion [W. Zurek, Physica Scripta, T76, 186 (1998)] that 
the chaotic tumbling of the satellite Hyperion would become non-classical within 20 years, but 
for the effects of environmental decoherence. The dynamics of quantum and classical probability 
distributions are compared for a satellite rotating perpendicular to its orbital plane, driven by the 
gravitational gradient. The model is studied with and without environmental decoherence. Without 
decoherence, the maximum quantum-classical (QC) differences in its average angular momentum 
scale as h 2 ^ 3 for chaotic states, and as h 2 for non-chaotic states, leading to negligible QC differences 
for a macroscopic object like Hyperion. The quantum probability distributions do not approach their 
classical limit smoothly, having an extremely fine oscillatory structure superimposed on the smooth 
classical background. For a macroscopic object, this oscillatory structure is too fine to be resolved 
by any realistic measurement. Either a small amount of smoothing (due to the finite resolution of 
the apparatus) or a very small amount of environmental decoherence is sumcient ensure the classical 
limit. Under decoherence, the QC differences in the probability distributions scale as (h 2 / D) 1 ^ 6 , 
where D is the momentum diffusion parameter. We conclude that decoherence is not essential to 
explain the classical behavior of macroscopic bodies. 

PACS numbers: 03.65.Sq, 03.65.Yz, 05.45.Mt 



I. INTRODUCTION 

Quantum mechanics (QM) is a more fundamental the- 
ory than is classical mechanics (CM). This fact does not 
diminish the utility of CM in describing the behavior 
of macroscopic objects. But theoretical consistency de- 
mands that the observed classical phenomena should also 
emerge from QM in an appropriate limit. However, a de- 
tailed understanding of this quantum-to-classical (QC) 
limit is remarkably difficult to achievc, and the origin of 
classical behavior is even subject to a degree of contro- 
versy. There has been some confusion as to which QM 
structures should be compared with classical predictions, 
and hence to what the appropriate criterion for classical- 
ity should be. The role that decoherence due to environ- 
mental perturbations might play in establishing the QC 
limit is also controversial. 

A criterion for classicality that has often been used 
is based on Ehrenfest's theorem 0. If the size of the 
quantum state is small compared to the scale on which 
the potential energy varies, then the centroid of the state 
will approximately follow a classical Newtonian trajec- 
tory. The time duration (and range of other relevant 
parameters) within which this condition holds is referred 
to as the Ehrenfest regime. Since any wave packet will 
eventually spread until it reaches the size of the system 
(harmònic oscillators being a unique exception) , it follows 
that this criterion for classicality will eventually fail, no 
matter how macroscopic the system may be. If the break 
time when this occurs were as long as the age of the so- 
lar system, there would be no cause for concern. But in 
chaotic systems, the size of a wave packet grows exponen- 
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tially, and the break time can be quite small. A striking 
example, given by Zurek 0, is the chaotic tumbling of 
Hyperion (a moon of Saturn) , for which the Ehrenfest cri- 
terion for classicality will fail in less than 20 years. Zurek 
argues that environmental decoherence can remove this 
paradox and restore classical behavior to Hyperion. This 
claim is certainly not correct, as long as the Ehrenfest 
criterion for classicality is used. However, we propose 
another criterion for classicality, within which the role of 
decoherence will be reassessed. 

The origin of the above paradox consists in a failing to 
take proper account of the statistical nature of QM. It is 
not correct to identify the trajectory of a body with the 
motion of the centroid of the wave function. QM does 
not describe the actual observed phenomcnon, but only 
the probabilities of the various possible phenomena. We 
should, therefore, compare quantum probabilities with 
classical probabilities. This leads us to define the Liou- 
ville regime of quantum- classical correspondence, within 
which the quantum probabilities are approximately equal 
to the classical probabilities that satisfy the Liouville 
equation. There is no requirement that the probabil- 
ity distributions be narrow, and so the Liouville regime 
of classicality is usually much larger than the Ehrenfest 
regime. 

The superiority of the Liouville criterion for classical- 
ity to the Ehrenfest criterion has been demonstrated in 
several ways. One way to see this is through the cor- 
rection terms to Ehrenfest's theorem. When the width 
of the state is small but not negligible, these corrections 
can be obtained as a series involving the variance and 
higher moments of the position probability distribution 
This series has exactly the same form for the classi- 
cal Liouville probability distribution. Hence Ehrenfest's 
theorem merely asserts that if the (quantum or classical) 
probability distribution is sufficiently narrow, its centroid 
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will follow a Newtonian trajectory. To leading order, the 
deviations from the Newtonian trajectory are the same in 
both the quantum and classical cases, and so those devia- 
tions are not primarily quantal in origin. The magnitude 
of the deviations from Ehrenfest's theorem depends pri- 
marily on the width of the initial state in configuration 
space, and has no systematic dependence on h |4j. How- 
ever, the (much smaller) differences between the quan- 
tum and classical probabilities scale as ti 2 , indicating that 
they are truly of quantal origin. 

Our interest in Hyperion was stimulated by Zurek's 
provocative paper ■ He begins with an estimate of the 
break time, Te, beyond which the Ehrcnfcst criterion of 
classicality will fail. Since the width of a wave packet in 
a chaotic system grows exponentially with the Lyapunov 
exponent A, the time taken for it to reach the scale L 
over which the potential varies (typically of order of the 
system size) will be Tg = À _1 ln(X/Axo), where Axq 
is the initial width of the wave packet. Zurek chooses 
Axq = Apo / fi-, with the width of the momentum distribu- 
tion Apo being estimated from thermal fluctuations, and 
thereby obtains the now-familiar result that the break 
time Te scales as ln(fi, _1 ). This Te can be quite short, 
even for a macroscopic system, and because of the log- 
arithmic dependence, it is not sensitive to detailed as- 
sumptions about the initial state. Thus, according to the 
Ehrenfest criterion for classicality, the tumbling motion 
of Hyperion should long ago have ceased to be classical. 

This paradoxical conclusion is not affected by includ- 
ing environmental decoherence. Decoherence converts a 
pure state into a mixed state, but it does not produce 
localization of the position probability density, and so it 
has no significant influence on the breakdown of Ehren- 
fest's theorem. (Indeed, the diffusive term that describes 
the effect of the environment in the master equation for 
the density matrix will have a slight delocalizing effect.) 
Therefore, if the Ehrenfest criterion were the sole cri- 
terion for classicality, we would still be faced with the 
paradoxical conclusion that macroscopic bodies like Hy- 
perion should be grossly non-classical. 

In @- EÜ| Zurek next considers the equation of mo- 
tion of the Wigner function, which has the form of the 
classical Liouville equation plus a series of fi-dependent 
terms, called the Moyal terms. The limit of the Liouville 
regime of classicality will presumably be reached when 
the Moyal terms have a significant effect. Unfortunately, 
Zurek does not distinguish between the Ehrenfest and Li- 
ouville regimes, and denotes both break times as tu. This 
is a serious confusion, since they are both conceptually 
and numcrically distinct. 

While the time limit Te of the Ehrenfest regime is 
easy to estimate reliably, an analogous limit for the Liou- 
ville regime is much more difficult to obtain. Zurek uses 
heuristic arguments to claim that decoherence tends to 
counter the effects of the Moyal terms, and he obtains a 
Liouville break time similar in form to Te 5}. We con- 
sider this conclusion to be doubtful for several reasons. 
Habib et al. [|| have shown that, in order for the non- 
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FIG. 1: Orbit of satellite spinning about the z-axis perpen- 
dicular to the orbital plane. 9 denotes the position of the 
satellite on the orbit, and <f> is the orentation of the satellite 
with respect to the semi-major axis of the orbit. 



negativity of the density matrix ( "rho-positivity" in their 
terminology) to be preserved, the Moyal terms must have 
subtle effects that are not counteracted by decoherence. 
Full numerical computations for a driven system Q and 
for an autonomous system Q have found that, although 
the differences between quantum and classical averages 
of observables have a brief period of exponential growth, 
these differences reach a saturation value, about which 
they fluctuate irregularly. This saturation value is much 
smaller than the system size (unlike the deviations from 
Ehrenfest's theorem), and it tends to zero as some small 
power of H. The break time estimated by Zurek would 
be relevant only if the exponential growth of these QC 
differences continued until they reached macroscopic size. 

In this paper we perform classical and quantum com- 
putations for the chaotic tumbling of an object like Hy- 
perion, which complement the published classical theory 
li- The actual value of the de Broglie wavelength is, of 
course, much too short to be treated in a numerical inte- 
gration of the Schròdinger equation. So we cast the equa- 
tion into dimensionless form, and solve it for a range of 
the dimensionless h parameter. These results lead to scal- 
ing relations, from which we can extrapolate to estimate 
the vàlues appropriate to Hyperion. We study the system 
with and without environmental decoherence, discuss the 
conditions needed to ensure classicality, as well as exam- 
ine whether there is a qualtiative difference between the 
classical limit of regular and chaotic initial states. 



II. MODEL 

Our model of Hyperion's rotation was first suggested 
in 1988 by Wisdom @. It assumes that Hyperion's cen- 
ter of mass travels in an elliptical orbit about Saturn, 
and that its orbit is unaffected by its rotation. However, 
since Hyperion is an extended object, Saturn's gravita- 
tional field is not constant over its volume. Since its mass 
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distribution is not spherical, the variation in the gravi- 
tational field can produce a net torque. To lowest order 
in a multipole expansion of the mass distribution, this 
torque depcnds on the quadrupole moment of the mass 
distribution, and for simplicity we neglect all higher order 
moments. 

It should be noted that this configuration is attitude 
unstable, and so small inclinations of I3 towards the plane 
of the orbit will tend to grow. However this simplifying 
assumption makes both the classical and the quantum 
mechanical computations feasible. 

The coordinate system is shown in Fig. ^ The space- 
fixed x-axis is along the semi-major axis of the orbital el- 
lipse, and the z-axis is perpendicular to the orbital plane. 
The angle 9 denotes the position of the satellite in the 
orbit. The axis of the smallest moment of inertia (li), 
makes an angle <j> with respect to the x-axis, hence the 
angle between the body axis of I\ and the radius vector 
f is (j) — 9. The largest moment of inertia 7 3 is parallcl 
to the z-axis. The canonical coordinates for this system 
are the angular momentum and the orientation of the 
satellite {L z ,(j)}. 

The coupling of the gravitational field to the satellite 
is obtained by a Taylor expansion of the potential about 
the satellite's center of mass, 



A. Classical Equation Of Motion 

It is convenient to express the equation of motion 
in terms of dimensionless variables. We introduce the 
anisotropy parameter, 



(h - h 



(5) 



a dimensionless time (in units of the orbital period), 

r = t/T (6) 

and a dimensionless angular momentum J z in terms of 
the dimensional angular momentum L z 



J, 



L 7 T 



(7) 



To estimate a for Hyperion, we use the observed 
lengths of its principle axes (410 ± 10, 260 ± 10, 220 
± 10 km) [TU, and assume that it is an ellipsoid of uni- 
form mass density. Hence 
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Here Xi refers to the distance along the i'th space-fixed 
axis from the center of mass of the satellite. 
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= diag 
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2Gm —Gm —Gm 



(2) 



Here m is the mass of the gravitational source (Saturn), 
and r is the distance from the source to the satellite. The 
first order term in Eq. vanishes because the expansion 
is about the center of mass, and the second order term is 
related to the moments of inertia tensor, 



p[x 2 Sij — XiXj]d x (3) 



Using Kepler's third law, which states GM/a 3 
4-7T 2 /T 2 ,' the Hamiltonian becomes 



Here T is the orbital period, a is the length of the semi- 
major axis of the orbit, L z is the angular momentum 
about the z-axis, I3 = I zz is the moment of inertia for 
rotations in the orbital plane, and r(t) and 9{t) are the 
orbital coordinates of the satellite, which are functions of 
the period T and the eccentricity e. These functions are 
found by numerically integrating the equations of motion 
for the center of mass, using the code provided in [Í"(Íj| . 



Here n is half the length of the i th principle axis of the 
ellipsoid. The other moments of inertia are obtained by 
cyclically permuting the indices. Substituting Eq. (jSJ) 
into Eq. © yields 



(9) 



Hence a — 0.43 ± 0.04. In this work we used a slightly 
larger value, a = 0.5, because it leads to a more purely 
chaotic motion, whereas for a = 0.43, e = 0.1 there are 
large regular islands embedded in the chaotic sea. We 
wish to compare chaotic motions with regular motions, 
and the differences would be obscured by a mixed phase 
space. 

Following Wisdom || , we obtain the equation of mo- 
tion (in dimensionless variables) to be 



r(r) 



asin(2[0 - 6>(t)]) (10) 



B. Quantum Mechanics 

The quantum mechanics will be solved by integrating 
the Schròdinger equation in angular momentum repre- 
sentation. The state vector is written as 



\m) =$> m (t)|m) 



(11) 
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with \m) being an angular momentum eigenstate. The 
matrix elements of the Hamiltonian are 



- h 2 m 2 c m ^ 37r 2 fa 



2U 



i r 2 * 

(J 2 -I x )— / ( e i(«-™+2)<A e -2i0 + ^(n-m^^iS^ 

(12) 

Using l|12|) and (|llfl . the matrix equation (to|/J|?/;) = 
ï?ï(to| J^|?/;) becomes 



,dc m (t) h 2 m 2 c m 3n 2 /a\ 3 
in — 



dt 



2I 3 



3* £)<*-*> 



(c m+2e 2ie W+c m _ 2 e-^^) (13) 
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value of ií must be chosen so that this range includes 
all of the vàlues of J z that have significant amplitudcs 
in the time-dependent state. By examining phase-space 
diagrams for the classical distributions, we found that 
| J z \ < 20 for all time, and so K = 20//3 was sufficient to 
contain the quantum distribution. 

The initial classical probability distributions are cho- 
sen so that they match the angular momentum and an- 
gular distributions for the initial quantum state. Because 
the initial state is a minimum uncertainty state with fixed 
width in angular momentum, its width in angle is pro- 
portional to (3. Thus (3 (dimensionless h) enters into the 
classical calculation to ensure that the initial quantum 
and classical states correspond to each other. 



III. RESULTS FOR A NON-CHAOTIC STATE 



In addition to the dimensionless parameters r and a, 
we now introduce a dimensionless h parameter, 



(14) 



The dimensionless Schròdinger equation then becomes 



.dc m 



j3m 2 c n 



3tH a 



2 (3 \r(r) 



{c m+2 e 2 ^) 



c m _ 2 e 



-2í6{t)\ 



(15) 



A peculiar feature of Eq. I|15[) is that the coefficient c m 
depends only on c m +2 and c m _2, therefore the even c m 
cannot interact with odd c rn . This coupling arises from 
the invariance of the Hamiltonian under rotations by n. 
But octapole and other odd moments are not invariant 
under rotations by n, so this symmetry is an artifact of 
the model. 



C. Initial State 

The initial quantum state is chosen to be a Gaussian 
in angular momentum, 



IV>> 



exp 



Qgm - Jq) 2 
2S 2 



(16) 



Here 0m is a dimensionless angular momentum, Jo is 
the average of the dimensionless angular momentum in 
the state, S is its Standard deviation, and 4>o is the central 
angle of the initial state. These parameters will be varied 
to ensure that the initial states are in regions of phase 
space that are either purely chaotic or purely regular. 

In principle the sum is from m = -oo to +oo, but 
in practice it is restricted to a range {—K-- -K}. The 
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FIG. 2: Poincare Section for e 
denotes a typical initial state. 



0, a = 0.5. Black circle 



The classical limit of the quantum tumbling of a satel- 
lite will now be examined for a non-chaotic state, to de- 
termine whether there is a qualitative difference between 
chaotic and non-chaotic systems in their approach to clas- 
sicality. 

Non-chaotic motion is ensured by choosing a circular 
orbit: e = 0, r(r) = a, 9(t) = 2t:t. The time dependence 
in Eq. Qlü[) can be transformed away by the substitution 
<I> — <f) — 2itt, yielding an integrable equation of motion, 



$ = -67r 2 asin(2$) 



(17) 



Fixed points for this equation oceur at the angles $ = 



— These fixed points describe motions in which 
Hyperion presents the same face to Saturn at all times. 
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FIG. 3: QC differences in (J z ) vs r, for several j3, with e = 0, 
a = 0.5. For /3 = {0.0125, 0.05, 0.5} the statistical errors are 
a m = {0.0002, 0.0007, 0.002}. 



FIG. 5: Scaling of the early time QC differences in (J z ) with 
(3, for e = 0, a = 0.5. Showing a f3 2 dependence. 



A. QC Differences in (J z ) 
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FIG. 4: QC differences in (J x ) for = 0.0125, with e = 0, 
a = 0.5, a m = 0.7 x 10~ 3 



The stable fixed points correspond to the smallest mo- 
ment of inertia pointing towards Saturn. 

The initial state was chosen to be far from the unstable 
fixed point. It is centered at Jo = 4, with a Standard 
deviation in J z of o~ = (see Fig. and a central 
anglei 4>q equal to zero. 



The classical probability distributions are found by 
time evolving a finite ensemble of systems, using Eq. 
(110(1 . The distributions of J z and 4> are found by ran- 
domly choosing the angular momentum and orientation 
of each membcr of the ensemble from probability distri- 
butions in J z and cj> that correspond to the initial quan- 
tum state. 

The finiteness of the ensemble leads to statistical er- 
rors, which may be reduced by increasing its size. The 
Standard deviation of the fluctuations in the mean is 

O-m = "4= (18) 
y/71 

Here n is the number of members in the ensemble, a is 
the Standard deviation of the distribution. Any difference 
between the computed mean vàlues of the quantum and 
the classical variables is not significant unless it is larger 
than cr m . Ensembles of 1,000,000 to 20,000,000 partides 
were used to ensure that the typical QC differences are 
greater than a m . 

As (3 — > the QC differences become smaller, and thus 
a larger ensemble is needed to reduce the statistical errors 
below that level. Hence different ensemble sizes were 
used for different vàlues of (3 in Fig. |3| The ensemble 
sizes were chosen so that a m = {0.0002, 0.0007, 0.002} 
for (3 = {0.0125,0.05,0.5}. 

Ensembles were evolved for several vàlues of (3, ranging 
from /3 = 0.5 to f3 = 0.002. For (3 < 0.01 the QC differ- 
ences were far smaller than <r m for any computationally 
feasible ensemble sizes, so no data will be presented for 
/3 < 0.01. 

A plot of the QC differences in (J z ) is shown in Fig. 
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FIG. 6: Variation of \qm — cl\i [Eq. I|19^ ] with time and fi for 
the nonchaotic state (e = 0, a = 0.5). Each classical ensemble 
has 1,000,000 members. 

13 In this and similar figures, any QC differences smaller 
than a m should be ignored, since they are dominated 
by statistical errors. The QC difFerences oscillate on the 
scale of the driving force, and only the envelope of these 
oscillations is of interest. From Fig. [21 it is apparent that 
at early times the envelope of the QC difFerences grows as 
t 2 . For longer times the envelope oF the QC difFerences 
is oscillatory, as can be seen in Fig. 0] Such recurrences 
are typical for non-chaotic systems [Tj|. Fig. [3] shows 
that, for fixed times, the QC difFerences in (J z ) scale as 
/3 2 . This result is similar to that found For some other 
systems Q. 

B. QC Differences in Distributions 

The differences in (J z ) alone are insuficient to fully de- 
scribe the differences bctwccn quantum and classical sys- 
tems because two different probability distributions can 
have the same mean but different variances and higher 
moments. We shall now examine the difFerences between 
probability distributions, and how they scale with /3. 

Since the angular momentum distributions are dis- 
crete, one can regard them as vectors, and measure the 
difference between the quantum and classical probability 
vectors by the 1-norm, defined as 

\qm-d\ 1 = ^\P cl {m)-P qm {m)\ (19) 

ÏTi 

This probability distribution is normalized so that 
~Y^ m P{jn) = 1. Alternatively, one can define a probabil- 
ity density, which is normalized so that J P(hm)d(hm) = 



1. Then the 1-norm of the probability densities takes the 
form 

/oo 
\P c i(hm) - Pg m (hm)\d(hm) (20) 
-oo 

These two forms are equivalent because P(hm) — 
P(hm)/h, and the additional Factor of H is cancelled by 
the Factor oF h in the integral. 

Fig. [ïïjshows that the QC differences in the probability 
distributions do not tend to zero as f3 — > 0. This lack of 
pointwise convergence of the quantum probability distri- 
butions to the classical Hmit has also been observed for 
other systems, such as a particle in a box and the kickcd 
rotor 13]. In these one-dimensional driven system, the 
quantum probability distributions develop a fractal-like 
structure, and only the smooth background converges to 
the classical probability distribution. We will show in 
section IV Cl that a similar result holds for Hyperion. 

IV. RESULTS FOR A CHAOTIC STATE 




FIG. 7: Poincare Section for a chaotic state, a = 0.5, e = 0.1. 
Black circle denotes a typical initial state, with J z = 10 and 
8 = 0.5 

In this section the rotation of a satellite is investigated 
for a chaotic state. The previous value of a — 0.5 is used, 
but now the eccentricity is taken to be Hyperion's value 
of e = 0.1. The computation is carried out as in section 

nu 

The initial state is centered at dimensionless angular 
momentum Jo = 10, with a Standard deviation of a = 
0.5, and a central angle 4>q — 0. This state is in the 
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FIG. 8: QC differences in (J z ) vs r for a chaotic state, a = 
0.5, e = 0.1. For = {1 x 10" 3 ,3 x 10~ 3 ,8 x 10" 3 } the 
statistical errors are a m = {1.4 x 10~ 4 , 1.4 x 10~ 4 , 2.2 x 10~ 4 }. 
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FIG. 9: QC differences in (J z ) vs (3, for a chaotic state before 
saturation is reached, showing a /3 2 dependence. 



chaotic sea, far away from any regular torii, as can be 
seen in Figure0 The maximum Lyapunov exponent for 
the chaotic sea is A = 0.85. 

This state was evolved for several periods of the driving 
force, and the differences between the quantum and clas- 
sical results were computed. In Fig. |S]the QC differences 
in ( J z ) are initially dominated by statistical errors, which 
are approximately a m = 2 x 10~ 4 . The QC differences 
in (J z ) grow exponentially with time until they saturate 
at t « 6. This saturation occurs when the classical tra- 



jectòries ergodically fill the chaotic sea. For r > 20 the 
classical ensemble saturates at (J z ) ps 8.2. The quantum 
value of (J z ) also saturates at approximately the same 
value, but with irregular fluctuations superimposed. This 
suggests that the QC differences here are dominated by 
quantum fluctuations, once the probability distributions 
have saturated the chaotic sea. 



A. QC Differences in (J z ) for Early Times 

To determine how the QC differences scale with /?, we 
varied (3 with r fixed at the times of the peaks in Fig. 
[S] From Fig. [ü]it can be seen that the QC differences in 
(J z ) scale as (3 2 ; the same sealing as was found for the 
non-chaotic states. This (3 2 sealing was previously found 
by a different method for other systems [4( . The classical 
ensemble sizes were chosen so that a m — {1.6x 10 -4 , 1.6x 
10~ 4 , 2.2 x 10~ 4 } for j3 = {lx 10~ 3 , 3 x 10~ 3 , 8 x 10~ 3 }. 

In the initial growth region of Fig. [SJthe QC differences 
in (J z ) vary with time as 



\{ j z)q ~ (Jz)c\ oc e 



2.9r 



(21) 



for r — 2 between 2 and 5.5. The exponent in Eq. I|21l) 

appears to be independent of the value of (3. The expo- 
nent is greater than 2A, implying that the QC differences 
grow at a rate that is greater than the classical Lyapunov 
exponent. Similar results have been obtained for some 
other systems 0, 0, Q • 

The exponential growth eventually saturates. This ces- 
sation of exponential growth of the QC differences in (J z ) 
is relevant to Zurek's argument ,2] that, absent deco- 
herence, the QC differences for Hyperion should reach 
macroscopic size within about 20 years. That argument 
implicitly assumes that the QC differences will continuo 
to grow exponentially until they reach the size of the 
system. However, we find that not to be the case. 



B. QC Differences in (J z ) for the Saturation 
Regime 

The maximum QC differences occurs in the saturation 
regime. If these differences converge to as (3 — > then 
the classical limit will be reached for all times, and there 
will be no break time beyond which QC correspondcncc 
fails. 

At the beginning of the saturation region (Fig. [S] and 
Fig. 115(1 . the QC differences in (J z ) reach a maximum, 
before decaying to a saturation level, about which they 
fluctuate irregularly. Because of this fluctuation in the 
saturation regime, we calculate the time average of the 
QC differences. Here \{J z )q — (J z )c\ was averaged over 
r from 20 to 100. As shown in Figure lTTl these averaged 
QC differences tend to scale as /3 2 ^ 3 . 

For sufficiently small (3, the peak QC differences also 
seales as /3 2 / 3 (see Fig. II Up . This sealing also was found 
for the maximum QC differences in a model of coupled 
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FIG. 10: Maximum QC differences in (J z ) versus j3 for a 
chaotic state. e = 0.1, a — 0.5. Suggests that this maximun 
difïerence scales as (5 2 ^ 3 . 
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FIG. 11: QC differences in (Jz), averaged over r from 20 
to 100. Initial state is in the chaotic sea. e = 0.1, a = 
0.5. Suggests that (\{Jz)q — (Jz)o\) oc (3 2 ^ 3 in the saturation 
regime. 



pcndulums [é| , so it might be genèric for chaotic systems 
in the saturation regime. 



FIG. 12: Quantum probability density for r = 40.0, (5 = 0.002 



probability distributions are shown in Figures ir^l and UÜl 
We use the quantity \qm — cl\i (defined in Eq. I|19|l ) as 
a measure of the QC differences in the probabilities. 

As can be seen from Fig. 1141 \qm — cl\i increases with 
time before saturating, but fails to converge to as f3 — > 
0. Since pointwise convergence does not occur, ncithcr 
for the chaotic nor for the non-chaotic states, this lack of 
pointwise convergence is not a result of chaos. 

Most of the QC differences in the probabilities occur 
on a very fine scale, and a modest amount of smoothing 
is sufhcient to cause the quantum probability distribu- 
tions to better approximate the classical results. Fig. 
EUshows that the differences between the two distribu- 
tions are dramatically reduced by smoothing them over 
a small width. This smoothing process is discussed and 
compared to environmental decoherence in section lV Cl 

To summarize the results of this section, for early times 
the QC differences in (J z ) scale as (3 2 , and increase ex- 
ponentially with time. The exponential growth ceases 
when the probability distributions saturate the chaotic 
sea. Both the maximum vàlues and the saturation levels 
of the QC differences were found to scale as /3 2//3 . A small 
amount of smoothing can dramatically reduce the QC 
differences in the probability distributions, since most of 
the differences come from very fine scale structures in the 
quantum probability distributions. 



ENVIRONMENTAL EFFECTS 



C. QC Differences in Probability Distributions 

The probability distributions contain much more in- 
formation than do the averages of observables. Thcsc 



Interaction with the environment leads to decoherence 
and dissipation. Decoherence is a quantum effect that 
causes interference patterns to decay. The time scale 
upon which this happens is model dependant, and for 
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FIG. 13: Quantum and classical probability densities for r = 
40.0, P — 0.002. Both quantum and classical densities are 
convolved with a triangular filter of width 0.25 in J z 



FIG. 15: Variation of (J z ) with r with and without the ran- 
dom potential with a/V c h = 0.024, r c = 0.01. (chaotic state, 
P = 0.05) 
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some systems there is no single decoherence timescale 
[Ï5| . Strunz et al. suggest that for the ràpid deco- 
herence expected in macroscopic bodies, the decay time 
varies as a small power of h, and is not sensitive to the 
system Hamiltonian. 

Dissipation is a classical effect which results in diffu- 
sive spreading of the probability distributions. Oftcn, 
the timescale for dissipation is much longer than the 
timescale for decoherence, and dissipation is insensitivc 
to h, unlike decoherence. But both effects are present 



together, and it is not always easy to separate them. 

The effect of the environment on a quantum system is 
often treated by a master equation that has non-unitary 
time evolution. Because an initially pure state can evolve 
into a mixed state, it is necessary to compute the density 
matrix, which requires much greater storage than does 
the computation of a state vector. The requirements for 
storage and computation time scale like K 2 , where K is 
the number of basis vectors needed to store a state vector. 
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FIG. 17: \qm — cl\i vs r for j3 = 0.05, for varying D = a 2 T c /6. 
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FIG. 18: |gm — ci|i vs r for varying /3, with r c = 0.01 and 
cr/V c h = 0.012, for a chaotic state. 



An alternative method is to perform n evolutions of Eq. 
(|15fl with a different realization of the random potential 
added for each run. Averaging the probability distribu- 
tions that result from the each of the n runs is physically 
equivalent to tracing over the environmental variables. 
The advantage of this method is that the computational 
resources for each run scale as K, rather than K 2 for the 
master equation. On the other hand, to achieve good 
accuracy, a large number of realizations of the random 
potential must be considered, in order to reduce the sta- 
tistical errors in the quantum calculation. However the 
number of realizations of the random potential that was 



needed to get sumcient accuracy was considerably less 
than K, so this method was much more computationally 
efficient than integrating the master equation. 

A stochastic potential is included, in both the quan- 
tum and classical mechanics, to model the effect of the 
environment on the satellite. The simplest stochastic po- 
tential that yields a random torque is 

H mt = V R(t) cos(0) (22) 

Here R(t) is a correlated random function of zero mean 
and unit variance, Vq is the amplitude of the random 
potential, and r c is its correlation time. A correlated 
random function is used because the fluctuations in the 
environment do not oceur instantly, but rather they oceur 
and decay on some time scale r c . The correlated random 
sequence R(t) can be constructed from an uncorrelated 
sequence, as is shown in Appendix A. The results are not 
sensitive to the exact form of Eq. 1)22(1 . and qualitatively 
similar results were obtained when cos(4>) was replaced 
by cos(2cf>). As is shown in Appendices B and C, the ef- 
fects of the environment are expected to depend mainly 
on the product ct 2 t c , rather than on the two parameters 
separately. Therefore we label the results by the momen- 
tum diffusion parameter, D = ct 2 t c /6, which is derived 
in Appendix B. 

The environmental perturbation should be much 
weaker than the tidal force on the satellite. Hence we 
compare the interaction potential Eq. (|22|l to the am- 
plitude of the tidal potential, made dimensionless by di- 
viding by h/T 2 , which is V ch = 3V2V 2 a, or V ch « 21 
for a — 0.5. In all cases reported in this paper, the 
environmental perturbation was so weak as to have no 
significant effect on the classical results, so its only sig- 
nificant effect is to produce decoherence in the quantum 
results. The same parameters as in the previous chaotic 
case were used, a = 0.5 and e = 0.1. 

Many realizations of the random potential were com- 
puted, and the results averaged, to get an aceurate mea- 
sure of the effects of the environment. We used 500 re- 
alizations to obtain results that are not strongly affected 
by statistical errors. Fig. H5lshows the QC differences in 
( J z ), with and without the random environmental poten- 
tial. The environment has no significant effect at early 
times, but in the saturation regime the QC differences are 
reduced. Since the primary effect of environmental deco- 
herence is to destroy fme-scale structures in the probabil- 
ity distributions, which do not affect averages like (J 2 ), 
this result may seem surprising. In fact, a typical trace 
of (J 2 ) versus time for a single realization of the ran- 
dom potential will look very much like that from a run 
without the random potential in Fig. 1151 But as time 
progresses, the oscillations in (J z ) for different realiza- 
tions of the random potential tend to get out of phase 
with each other, and the decreased amplitudes of the QC 
differences in Fig. 1151 are due to the averaging over the 
many different realizations of the random potential. 

In Fig. ^jl 100 realizations of the interaction potential 
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FIG. 19: Quantum and Classical probability distributions at r = 40, with (5 — 0.0125, o /V c h = 0.012, and r c = 0.01 for the 
chaotic state. In (b) the sòlid lines denote the results with the environment, and squares without the environment, showing 
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FIG. 20: Decay times Td of \qm — cl\\ vs the correlation time r c 
of the perturbing environment, with (3 — 0.05, <r/V c h = 0.012, 
in the saturation regime (a = 0.5, e = 0.1) 



were sufíicient to find the maximum QC diffcrcnces in 
\qm — cl\i (Eq. However, 700 realizations of the 

interaction potential were insufficient to resolve the QC 
differences in \qm — cl\\ in the saturation regime, and so 
for computational reasons, these will be estimated rather 
than directly computed. The variation of these differ- 
ences with (3 and the environmental parameters can be 
seen in Figures. 1171 and ITH1 

It can be seen in Fig. ^5]that, with the inclusion of 
the environmental perturbation, the quantum probability 



distribution is much closer to the classical distribution 
than without the environment (compare Fig. I12f> . 



A. Environmental Effects in the Saturation Regime 

In the saturation regime, the QC differences in \qm — 
cl\i (Eq. lfTÜ|) ) appear to decay exponentially from their 
maximum value to a saturation level (see Fig. 1171 and 
I18|l . The rate of this decay is a duffusion time, and is not 
the decoherence time, as will be shown shortly. From 
Fig. [20| it is clear that, for sufficicntly largc vàlues of r c , 
the decay time no longer depends on r c , but settles at 
Td ~ 5.6. For sufficiently large a and sufficiently small 
(3, the decay time was found to also have approximately 
this same limit. 

To test whether this decay rate is governed by quantum 
mechanics, we compared two classical ensembles with dif- 
ferent initial vàlues of Jo (Jo = 10 and J = 11), and 
computed the 1-norm of the difference between them as 
a function of time. Fig. |2] shows that these initially 
different classical ensembles converge at a rate given by 
Td = 5.6. So, apparently, this time scale measures how 
quickly the differences between two different distributions 
decrease as they both grow to fill the chaotic sea. 

It is not clear from Fig. El an d El whether the QC 
differences in the probability distributions eventually de- 
crease to zero or reach a non-zero long-time limit. In Fig. 
1221 the long-time saturation level of \qrn — cl\i is plotted 
as a function of the number n of realizations of the ran- 
dom potential. In the limit n — > oo the QC differences 
approach a small value that appears to be slightly pos- 
itive. However that extrapolated limit is substantially 
smaller than the typical statistical errors for ensemble 
sizes of 1,000,000 to 2,000,000, and so is not significantly 
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FIG. 21: 1-norm of the difFerence between two classical en- 
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FIG. 23: Maximum vàlues of \qm — cl\i vs f3 2 /D. The points 
labeled /3, a, andr c represent data sets where /3, a, and r c were 
varied with the other two parameters held constant. 



different from zero. 
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with arguments presented by Pattanayak et al. [T^ , the 
data was found to collapse onto a single curve parame- 
terized by £ = (3 2 /D (Fig. EU). From a least squares fit, 
the scaling relationship was found to be 



max(|çm — cl\i) oc 



(23) 



This scaling as fi- 1 / 3 was also found for a coupled rotor 
model without decoherence Q . This result suggests that 
the /3 1 / 3 scaling found here might be genèric for systems 
with more than one degree of freedom, and also suggests 
that the pointwise convergence of the quantum proba- 
bility distribution to the classical distribution may occur 
because of other interacting degrees of freedom (not nec- 
essarily an external environmcnt). 



C. Effects of Decoherence vs Smoothing 



FIG. 22: \qm — cl\\ vs 7i~ 1//2 at r = 40. n is the number of 
realizations of the random potential. (a/V c h = 0.012, r c = 
0.01 for the chaotic state) 



B. Scaling of the Maximum QC Differences 

The maximum value of \qm — cl\± for the quantum and 
classical probability distributions must depend on the 
three parameters /?, a, and r c . However, in agreement 



The effect of environmental decoherence is analogous 
to applying a smoothing filter to the quantum distribu- 
tions, but the two effects are not exactly the same. A 
small amount of smoothing does not affect the average 
vàlues of an observable, whereas decoherence can have 
an effect by dephasing the quantum fluctuation, which 
are then averaged over (see Fig. I15fl . This effect may be 
peculiar to systems with only one degree of freedom. 

Anothcr diffcrcnce is in the scaling of the QC differ- 
ences with the parameters. Figure |2*H shows \qm — cl\i 
in the saturation regime as a function of (3 and A s , 
where A s is the half width of the triangular filter used 
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to smooth the quantum distribution. In the saturation 
regime, \qm — cZ]i was found to depend on the single 
parameter /3/A s , where it obeys the scaling relation 



\qm-cl\t = 0.58 



0.44 



(24) 



Eq. I|24|l shows a scaling relationship similar to Eq. 
B23J1 . This may be understood as follows. Environmental 
perturbations cause momentum difFusion. This effect is 
proportional to v^D, and so A s should be compared to 
D 2 . Hence, in the late time limit, smoothing and deco- 
herence rely on a similar composite parameter. However 
the power laws are different in the two cases, and so the 
two processes are not entirely equivalent. 

In summary, the environmental perturbations were 
found to drasticly reduce the fine-scale structure in the 
quantum distributions. The measure of the QC differ- 
ences, \qm — cl\\ , was found to initially increase with r 
in a form similar to the results in Sec. IIVI After reaching 
a maximum value, \qm— cl\\ then decreased exponentially 
with time. The decay time was found to be a classical 
diffusion time, and not a decoherence time. The maxi- 
mum QC differcnces were found to scale as ([3 2 /D) 1 / 6 , 
where D is the momentum diffusion parameter. 



VI. CLASSICAL LIMIT FOR HYPERION 

A. QC Differences Without Environment 

Having calculated the QC differences for the chaotic 
rotation of a tidally driven satellite, and determincd how 
they scale with the relevant parameters, we shall now use 



this information to estimate the magnitude of quantum 
effects on Hyperion. In particular, we shall assess Zurek's 
claim 0, Uf that environmental decoherence is needed 
to ensure its classical behavior. We first examine the 
magnitude of the QC differences for Hyperion if the effect 
of the environment is ignored. 

First we must determine the dimensionless parameter 
(3 = Using Hyperion's mean density of p = 1.4 g 
cm -3 , and treating it as an ellipsoid with moments of 
inertia I3 > I2 > li, then 7 3 = 2.1 x 10 29 kg m 2 . Using 
the value of h = 1.05 x 10~ 34 J s, and the orbital period 
T = 1.8 x 10 6 s, then yields 



= 9.3 x 10 



-58 



(25) 



In Sec. IV B it was found that the maximum QC 
differences in (J z ) scale as /3 2 / 3 . Hence the maximum 
QC difference in the dimensionless angular momentum 
(J z ) for Hyperion should be approximately 5 x 10 -37 . 
So there should be no observable difference between the 
quantum and classical averages of angular momentum for 
Hyperion. 

This result contradicts Zurek's claim that, if decoher- 
ence is ignored, there should be a break time of no more 
than 20 years, beyond which the QC differences would 
become macroscopic. As was pointed out in the Intro- 
duction, if that break time is interpreted as the limit of 
the Ehrenfest regime, then it does not mark the end of 
the classical domain. But in [5j a break time of a simi- 
lar order of magnitude was estimated for the end of the 
Liouville regime. Both of those estimates were based on 
an exponential growth of the QC differences that occur 
in a chaotic state. Now the deviations from Ehrenfest's 
theorem do, indeed, grow exponentially until they reach 
the size of the system, as is needed for Zurek's argument 
to succeed. But, as was shown in Sec. IV(A), the expo- 
nential growth of the differences between quantum state 
averages and classical ensemble averages will saturate be- 
fore those differences reach the size of the system, and the 
saturation value seales with a small power of h. Hence, 
for the actual (small) value of H, the QC differences in the 
Liouville regime can remain small for all time, and there 
is no effective break time for the regime of classicality. 

The differences in (J z ) become vanishingly small in 
the classical limit, but this does not imply that the full 
quantum probability distribution converges to the classi- 
cal limit. We know that the quantum probability distri- 
bution will not converge to the classical distribution in 
a pointwise fashion. But we can ask what resolution is 
needed for a detector to be able to discriminate between 
these two distributions. Let us suppose that two proba- 
bility distributions are practically indistinguishable when 
\qm — cl\i < 0.01. Using scaling result in Fig. [31 we ffnd 
that a resolution A s of 1 part in 10 -60 rad/s is needed 
to resolve the two probability distributions. This sug- 
gests that it would be practically impossible to observe 
the quantum effects in the probability distributions, even 
without invoking environmental decoherence. 



14 



B. Environmental Effects on Hyperion 



VII. CONCLUSION 



There are many environmental perturbations that can 
affect the satellite: random motion of the partides within 
the satellite, random col·lisions with interplanetary dust, 
and random light fluctuations from the sun, to name a 
few. We shall consider the random collisions with dust 
partides as an example. 

To do this we treat the interplanetary dust as a dilute 
gas, and Hyperion as a sphere rotating about a fixed axis 
under the influence of random motion of the fluid. The 
dimensional momentum diffusion parameter is |18| 



D = lÇmkTR\r) 



(26) 



Using Eq. l)B4[l . the dimensionless momentum diffusion 
parameter D is 



D 



8iTkTRlr)T 3 



T 2 

i 3 



(27) 



Here T is the temperature, k is Boltzmann's constant, Rh 
is the radius of Hyperion, and rj is the kinetic viscosity of 
the dust fluid, which, following [n| , is calculated to be 



nmvL 
3a/2 



(28) 



Here v is the rms velocity of the dust partides, m is their 
mass, n is their number density, and L = l/nirr 2 is their 
mean-free-path, and r is the radius of a dust particle. 

The properties of interplanetary dust were measured 
by the Voyager space probes. The average number den- 
sity of partides near Saturn is n = 4 x 1CP 8 m~ 3 pfij . 
The average mass of the dust grains is estimated to be 
m = 10 -10 g, and their radius is about r = 10 m. The 
temperature in the vicinity of Saturn is about T « 135 

k plg. 

Using Eq. (J2BJ, treating Hyperion as a sphere of ra- 
dius Rh = 150 km, and using rj = 1.8 x 1CP 6 Pa s for the 
kinetic viscosity, we estimate the dimensionless angular 
momentum diffusion parameter to be D « 6.4 x 10 -50 . 
Even such a small value is sufficient to reduce \qm — cl\\ 
substantially. Using Eq. i|23|) , the order of magnitude of 
\qm — cl\i for Hyperion is found to be 10 -10 . This implics 
that the classical and quantum probability distributions 
will agree almost exactly for a large body such as Hyper- 
ion. Without the influence of the environment, the value 
of \qm — cl\i due to the very fine-scale differences between 
the quantum and classical probability distributions might 
be of order unity. But, of course, these differences would 
be impossible to resolve because they exist on such a very 
fine seale. So the effect of decoherence is to destroy a fine 
structure that would be unobservable anyhow. 



In this paper the regular and chaotic dynamics of 
a satellite rotating under the influence of tidal forces 
was examined, with application to the motion of Hype- 
rion. Quantum and classical mechanics were compared 
for both types of initial state, and the sealing with h of 
the quantum-classical (QC) differences was determined. 
The effect of the environment was modeled, and its ef- 
fect on the QC differences was estimated, so as to de- 
termine whether environmental decoherence is needed to 
account for the classical behavior of a macroscopic ob- 
ject like Hyperion. Two measures of the differences be- 
tween quantum and classical mechanics were examined: 
the QC difference in the average angular momentum, 
\(Jz)q ~ (Jz)c\, and the differences between the prob- 
ability distributions, \qm — cl\i (Eq. lfTÜ|) ). 

For early times, the QC differences in (J z ) grow in 
time as r 2 for the non-chaotic state, and as e 2 9T for the 
chaotic state. At longer times, the QC differences satu- 
rate for the chaotic state, but oscillate quasi-periodically 
for the non-chaotic state. The magnitude of the QC dif- 
ferences scale as /3 2 (dimensionless K) at early times, for 
both the chaotic and non-chaotic states. This 1 sealing 
persists for all times for the non-chaotic state. But the 
QC differences that oceurs in the saturation regime of 
the chaotic state scale as /3 2 / 3 . A similar sealing has also 
been observed for a model of two coupled rotors y|, so 
this result is not peculiar to the particular model studied 
in this paper. 

The value of the dimensionless h for Hyperion is (3 — 
9.3 x 10~ 54 , for which the /? 2 / 3 sealing relation prediets 
a maximum value for the QC difference in (J z ) to be 
5 x 10 -37 . Therefore, there is no need to invoke environ- 
mental decoherence to explain the classical behavior of 
( J z ) for a macroscopic object like Hyperion. 

Although the differences between the quantum and 
classical averages of observables becomes very small in 
the macroscopic limit, this need not be true for the differ- 
ences between quantum and classical probability distri- 
butions. Indeed, the quantum probability distributions 
do not converge pointwise to the classical probability 
distributions, for either the non-chaotic or the chaotic 
states. A modest amount of smoothing of the quantum 
distribution reveals that it is made up of an extremely 
fine-scale oscillation superimposed upon a smooth back- 
ground, and its is that smooth background that con- 
verges to the classical distribution. Similar behavior has 
been found for other one-dimensional systems |13|. This 
smoothing can be regarded as an inevitable consequcncc 
of the finite resolving power of the measuring apparatus. 
Alternatively, it may be impossible to observe the fine 
structure because of environmental decoherence. At the 
macroscopic scale of Hyperion, the primary effect of de- 
coherence is to destroy a fine structure that is anyhow 
much finer than could ever be resolved by measurement. 

When the environment was included, the results were 
found to follow a sealing relationship proposed by [Ï7| : 
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the maximum distance between the classical and the 
quantum probability distributions is proportional to 
(/3 2 / D) 1 / 6 . Here D is the momentum diffusion parame- 
ter (see Aupendix IBI) . This suggests that the quantum 
probability distributions will approach the classical dis- 
tributions pointwise as /3 — > 0, provided that D is non- 
zero. With environmental perturbations included, the 
QC differences in the probability distributions scaled as 
\qm — cl\i oc /3 1 / 3 . A similar scaling was also found for 
two autonomous coupled rotors This suggests that 
pointwise convergence of the quantum probability distri- 
bution to the classical valuc may bc typical for systems 
with more than one degree of freedom, and the lack of 
such convergence for systems with only one degree of free- 
dom may be pathological. The role of the environment, 
in the model of this paper, is then to cure this pathology 
by supplying more degrees of freedom. 

Taking D to be the momentum diffusion parameter 
for rotation of Hyperion due to collisions with the inter- 
planetary dust around Saturn, we find (estimated from 
Eq. l[HJ) that the maximum of \qm — cl\\ that Hyperion 
should exhibit should be of order 10 -10 . Thus decoher- 
ence would cause the quantum probability distribution 
to converge to the classical distribution in essentially a 
pointwise fashion. 

Coarse graining (due to the finite resolution power of 
the measurement apparatus) will also decrease the QC 
differences in the probability distributions. In the satu- 
ration regime, the measure \qm — cl\\ of that difference 
was found to be proportional to ((3/ A S ) 0A4 , where A s is 
the width of the smoothing filter. This shows that de- 
coherence and smoothing have similar effects. But they 
are not exactly equivalent, since their effects scale with 
somewhat different vàlues of (3. 

In conclusion, wc find that, for all practical purposes, 
the quantum thcory of the chaotic tumbling motion of 
Hyperion will agree with the classical theory, even with- 
out taking account of the effect of the environment. De- 
coherence aids in reducing the quantum-classical differ- 
ences, but it is not correct to assert that environmental 
decoherence is the root cause of the appearance of the 
classical world. 



APPENDIX A: CORRELATED RANDOM 
NUMBER GENERATION 

To describe the effect of environmental perturbations, 
we require a sequence of correlated random numbers. 
Generators for uncorrelated random variates are com- 
monly available, but algorithms for gcnerating a corre- 
lated sequence are not common. We show here how to 
generate a random sequence having a controlled amount 
of correlation from a Standard sequence of independently 
distributed random numbers. Let {r^} be such a se- 



quence, with zero mean and unit variance. 

(n) = (Al) 
(riTj) = % (A2) 

To generate a correlated sequence {Ri} from the un- 
correlated sequence, we simply form linear combinations, 

Ri+í = cR t + (1 - c)r l+1 (A3) 

where c is a chosen positive constant (c < 1), and R\ = 
ri . It follows from Eq. (|A3|) that 

i-l 

Ri+i = à-\ x + c m (l - c)r i+1 _ m (A4) 

m=0 

From this result, we can calculate the degree of correla- 
tion in our new sequence. Taking i > j, and using Eq. 
(|X4|) and Eq. (|A^|) . we obtain 

i-2 ]-2 

(RiRj) = J2J2 c m+m '(l - c) 2 (r l+1 _ m r J+1 _ m ,) + 

m=0 m'=0 

i-2 
m—O 

]T c m '(l - cjc^fa+x-^n) + c^- 2 {rf) 

m'=0 

(A5) 

Eq. (|A5I) can be simplihed using Eq. I|A2|I . Performing 
the resulting geomètric sums then yields 

{RiRi) = (^) à- 1 ^ c^) (A6) 
For j > 1 and c < 1, Eq. I|A6(1 becomes: 

(RiRi) « (^) (A7) 

This discrete sequence must now be converted into a 
function of time. Each Ri refers to the correlated random 
function at time t j . Taking the time interval between the 
random numbers to be Aí, then it is appropriate to de- 
fine a correlation time r c = Aí/| ln(c)| for the correlated 
random function, for which we have 

(R(0)R(r)) » (j^ e~^° (A8) 
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APPENDIX B: MOMENTUM DIFFUSION 
PARAMETER 

The momentum diffusion parameter (D) is needed to 
calculate the efïect of the environment on the system 0, 
Il6l l2ll | . In particular, the form of D is needed to show 
that the scaling result in apphes also to our model. 

Consider the random potential of the form V — 
Vq cos(0). The random torque is thcn F = Vo-R(í) sin(0). 
Here R(t) is a correlated random function, as defined in 
Appendix A. From this, we can find the momentum dif- 
fusion parameter through the relation 



s(t) 2 ) 



D = lim 

í^O t 



(Bl) 



The integral of the torque over time yield the angular 
momentum, hence the variance of the angular momentum 
under this random torque is given by 

{s 2 {t)) = % f f dt> dt" (R(t')R(t")) 



r 2 
1 



dt'dt" fl^ ) c i*'-*" < 

2 Jo Jo \l + c, 



\,- I ' C Ih2 > > <2 



1 + c 



{tj + t c t + t 2 c exp(-í/t c )) 



(B2) 



The quantity s(t) is the Standard deviation of the angu- 
lar momentum for a random walk under the influence of 
Eq. (E3), and (7 is the Standard deviation of the random 
potential. For t 3> t c we have 



(s 2 (t))^tt c V 2 



1 - c 
1 + c 



(B3) 



Using IBÏI and choosing the value c = 1/2, we obtain 
momentum diffusion constant D to be 



D = 



vtt 



L c 



(B4) 



In the body of this paper, we use a dimcnsionless mo- 
mentum diffusion parameter D. The relation between 
these two quantities is 



D = D 



rp3 

2J| 



become a function of the single parameter £, where 
£ = h a X b D c , for some powers a, b and c. Here D is 
the diffusion parameter, defined in Appendix B, and A is 
the classical Lyapunov exponent. 

The argument is similar to one presented in 0- It as- 
sumes that the differences between quantum and classical 
mechanics are due to the Moyal terms in the equation of 
motion for the Wigner function. These terms have the 
form 



\ ^ 



^ 2 2n (2n + 1) 

71—1 K ' 



■dl n+1 V{4>)d 2 p n+1 p w 



(Cl) 



where p w is the Wigner function for the state. For a 
chaotic system, the phase space distribution will develop 
very hne structures as it fills the accessible phase space, 
with the rate at which these fine structures develop be- 
ing governed by the Lyapunov exponent A. Since these 
terms depend on dp n+1 p w , they will become larger as 
time progresses and the fine structure grows. However 
the inclusion of environmental perturbations on system 
causes diffusion, which will limit the growth of the fine 
structure. When these effects balance each other, the fine 
structure is expected to have an equilibrium scale 0, ílTj 
given by 



(C2) 



If this equilibrium momentum scale is sufficiently large 
compared to h, then only the first order Moyal term 
should be significant. Under this assumption, the QC 
differences should be a function of this Moyal term, given 
by 




2 x 5/2 n -(3/2) 



£ = h z X >lz D 



(C3) 



p 



It is argued that the characteristic scale on which 
varies is S<p oc y/D/X. The characteristic variation of 
d^V(cj)) can be found by a Taylor expansion of d^V((f>) oc 
sin(0) about an arbitrary point 4>o — 0. Averaging the 
result over <j>o = then yields 



(C4) 




( B5 ) Inserting Eq. (Ü4j into Eq. (ÜSJl yields 



where r c is defined to be r c = t c /T, and a = VoT 2 /I 3 . 



APPENDIX C: DERIVATION OF SCALING 
PARAMETER 

Pattanayak et al. 01 suggested that at long times 
the QC differences in the probability distribution should 



(C5) 



Thus, once the growth fine scale structure of the proba- 
bility distribution reaches equilibrium with diffusion, the 
QC differences in the probability distribution should be 
a function of fi, 2 !? -1 . The results in Sec. Ivl confirm this 
conclusion. 
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